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A microscopic theory for rare-earth ferromagnetic hexaborides, such as Eui_a;Caa;B6, is proposed 
on the basis of the double-exchange Hamiltonian. In these systems, the reduced carrier concen- 
trations place the Fermi level near the mobility edge, introduced in the spectral density by the 
disordered spin background. We show that the transport properties such as Hall effect, magne- 
toresitance, frequency dependent conductivity, and DC resistivity can be quantitatively described 
within the model. We also make specific predictions for the behavior of the Curie temperature, Tc, 
as a function of the plasma frequency, uip. 



PACS numbers: 71.23.An, 75.47.Gk, 75.47.-m 

The series of compounds Hi-x-^x^e, where A is an 
alkaUne-earth metal such as Ca or Sr, and R a rare- 
earth magnetic ion, has recently attracted considerable 
interest, following a series of experiments which unveiled 
intriguing connections between their magnetic, transport 
and opticalproperties. These are cubic compounds where 
a divalent [ij lanthanoid occupies the central position on 
a cube, surrounded by eight Bg octahedra at each ver- 
tex. Boron atoms make up a rigid cage, held together by 
covalent bonds between neighboring B atoms. Our fo- 
cus will be on the results known for Eui_a;Caa;B6. EuBe 
is a ferromagnetic metal, ordering at Tc ~ 15 K, and 
characterized by a quite small effective carrier density jOf 
order of 10"^ per unit cell, at high temperatures 
Magnetism is found to arise from the half-filled 4/ shell 
of Eu, whose localized electrons account for the measured 
magnetic moment of 7 /ie per formula unit . The 

FM transition temperature is reported to decrease with 
increasing Ca content 0, Q and the totally substituted 
compound CaBg exhibits no magnetism. 

In this paper we propose a simple model that describes 
quantitatively the properties revealed by the experiments 
done in EuBg: (i) a precipitous drop in the DC resistivity 
just below Tc, with a change by a factor as high as 50 
between Tc and the lowest temperatures 0, 0| ; (ii) the 
large negative magnetoresistance observed near Tc; (iii) 
an increase in the number of carriers, by a factor of 2 - 3, 
upon entering the ordered phase, as evidenced by Hall 
effect 01; (iv) a large blue shift of the plasma edge, seen 
also for T < Tq, both in reflectivity, R{lo), and polar Kerr 
rotation (2, tl la| ! (v) a scaling of the plasma frequency 
with the magnetization. The above Hsted features consti- 
tute a case for the strong coupling of the magnetization 
to the transport properties. 

The effects of chemical doping with non-magnetic Ca 
are also considered, and the theory explains qualitatively 
the following experimental flndings: (1) with doping, x, 
the metallic regime, found in EuBg [x = 0), evolves to a 
semiconducting behavior above Tc [4, _6||; (2) just below 
Tc the carrier density increases by at least two orders 



of magnitude 01; (3) the plasma edge is visibly smeared 
while the corresponding resonance in the polar Kerr ro- 
tation is greatly attenuated 0,0|; (4) p{T,H) and ujp 
display an exponential dependence in the magnetization 
laUj; (5) there remains a signiflcant and rapid decrease 
of p{T) just below Tc, albeit by a smaller factor than in 
the undoped case. Our model predicts that the square of 
the plasma frequency, uj^, scales linearly with the Curie 
temperature, Tc, as doping x is varied. Furthermore, we 
also present a phase diagram T versus x which contains 
two phase transition lines: paramagnetic to ferromag- 
netic and metal to insulator. 

Theoretical understanding of the microscopic mecha- 
nisms responsible for the transport data is still contro- 
versial. EarHer local density approximation (LDA) cal- 
culations predicted semi-metallic character with an over- 
lap of conduction and valence bands at the Fermi level 
on the X point of the Brillouin zone (l^ IT^ . This 
appears to be consistent with the Shubnikov-de Haas 
and de Haas-van Alphen measurements 0, 0|. Fur- 
ther density functional theory (DFT) calculations showed 
the existence of a sizeable gap of ~ 1 eV. These results 
are corroborated by angle resolv ed p hotoemission spec- 
troscopy (ARPES) for magnetic [l^ and non-magnetic 
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hexaborides [Tg]. Some existing 
sume the semimetal scenario 
seems to rule out. 

We present a microscopic theory for EuBg, involving 
a single conduction band, consistent with the observed 
properties of these compounds, and make specific pre- 
dictions regarding the phase diagram for Eui_a;Ca2;B6. 
The presence of defects in the crystalline environment 
l2(l |2l| is believed to create a small carrier concen- 
tration in the conduction band. These carriers interact 
with the local Eu magnetic moments in a cubic lattice, 
via an s-f exchange coupling, leading to the well known 
Kondo lattice Hamiltonian: 
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where Si represents the 4/ Eu spin {S = 7/2) opera- 
tor at the site i and = ^/2J2a i3'4,a''^a,i3Ci^p is the 
the electron spin operator. Given the large value of 
S, the core spins are treated as classical variables pa- 
rameterized by the polar and azimuthal angles as Si = 
S {sin 0i cos (fii, sin 0i sine/) i, cos Oi). Based on the current 
literature, reasonable ranges for the model parameters 
are Uj = t ~ 0.1-1 eV and JS ~ 0.35-0.7 eV. 

The carrier concentration in EuBe {n ~ 10"'^) is very 
small, and therefore the mean kinetic energy is much 
smaller than the magnetic one, {K) ^ tn ^ JS. In 
this regime, the low energy physics of Hamiltonian (QJ 
can be obtained by taking the limit of J ^ oo [23|. The 
effective Hamiltonian is obtained by projecting out the 
states with the electronic spin parallel (for J > 0) or anti- 
parallel (for J < 0) to the local core spin S^. In either 
case, the result is nothing more than the double-exchange 
(DE) Hamiltonian (il: 

n=-Y,{\,c)c,+h. c) . (2) 

The effective hoping amplitude is ti^j = 
t cos (6*, /2) cos [9 J /2) -I- sin [6^ /2) sin {6^ /2) e~^('^'-'^^-) . 

The Hamiltonian Q describes a very interesting in- 
terplay between the magnetic and electronic degrees of 
freedom. The spin texture acts as a static disordered 
background for the electronic motion The ground 

state of the problem is obtained by minimizing the kinetic 
energy with the alignment of all the spins in the system. 
Thus, the electronic density of state (DOS), N{E,M), 
and the Fermi energy, Ep{M), change with magnetiza- 
tion, M, even if the number of electrons is constant. In 
this ferromagnetic state the transport properties are tied 
to the magnetic ones. Similar physics can be found in the 
context of the colossal magnetoresistance (OMR) man- 
ganites [i^. This is indeed an intrinsically disordered 
problem where the strength of the non-diagonal disorder 
is determined by M. A mobility edge, Ec{M), appears 
in the spectral density and strongly depends on M. 
In the absence of structural disorder, when the system is 
fully magnetized (M — 1) the mobility edge moves out 
of the band. One immediately realizes that in the para- 
magnetic phase (M = 0), if Ec{0) and Ef{M = 0) are 
comparable, the onset of ferromagnetism at Tc and the 
concomitant displacement of Ec towards the bottom of 
the band, have a major effect in the transport properties, 
by allowing more states to become delocalized and con- 
sequently increasing the conductivity. This was actually 
a concept that lingered for some time in the context of 
the manganites where the metal-insulator transition and 
CMR were assumed to stem from this effect In that 
case spin disorder alone cannot account for the experi- 
mental evidence since Ec{0) encloses less than 0.5 % of 
the states In the case of EuBg, however, the small 
number of localized states is comparable to the concen- 



1x10"' 




Temperature (K) 



Figure 1: Evolution of ne{T) assuming = 0.003 and 
Ew = 0.1 1. The inset shows the mobility edge Ec (circles) 
and mobility gap A (squares) as functions of the normalized 
magnetization M. 

tration of mobile carriers and thus this effect should have 
visible consequences. 

"Pure" EuBq. In order to investigate the aforemen- 
tioned, we calculate N{E, M) and Ec{M) for the Hamil- 
tonian l(2l for static spin configurations with specific val- 
ues of M. The full, self-consistent treatment of the 
magnetism of Hamiltonian Q is computationally too 
demanding for the system sizes required to study lo- 
calization effects. In this work we circumvent these 
difficulties by assuming uncorrelated spins in an exter- 
nal field and extracting the temperature dependence of 
the zero-field magnetization, M(T), from the experi- 
mental data 0. The recursion method |2^ was used 
to calculate N(E,M), and Ec(M) was located using 
the transfer matrix method {3flL 13 1| . Non-magnetic de- 
fects were accounted by a small magnetizationjndepen- 
dent shift, Ew, of Ec, thus giving effectively Ec{M) = 
Ey/ + Ec{M). The important quantity is the mobil- 
ity gap, A(M) = Ec{M) ~ Ef{M) < 0, the variation 
of which, at the transition, completely determines the 
number of extended (metallic) carriers, ne(M). Tuning 
Ef{Q) so that Tig = ne(0) — 0.003, as reported in the 
transport measurements above Tc, the variation of Ue is 
presented in Fig.Q] A(Af) is found to be almost linear in 
M: A(M) = Ao(l - aM) (see the inset in Fig.[TJ. The 
range of variation of rie in Fig. compares well with the 
data obtained from the Hall effect by Paschen et al. Q|. 

Another experimentally accessible quantity is the 
plasma frequency, ujp. Applying the Kubo formula for 
the optical conductivity, o-(w), to Hamiltonian (|2| the 
sum rule for its real part, <y' [t-o), can be extracted and 
combined with the optical sum rule [s^ leading to the 
plasma frequency for the model uf^ = — 47re^a^ {Ti) /3. 

In the reported experimental results Q , is obtained 
from the Drude contribution to the optical response of 
the material. Therefore, in order to compare with the 
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Figure 2: Result for Ljp{M) obtained from the sum rule (cir- 
cles) and experimental data (diamonds). The experimental 
{H — 0) points were obtained combining Wp{T) at zero field 
from ref. with the zero field M{T) from ref. Q. 



model predictions, we included only the contribution of 
metallic states in the calculation of the kinetic energy, 



a;2(M) ^ !Z(M)N{E,M)f{E)EdE. 



lation ijjp 
indeed verified 



The outcome of 
3| The linear ra- 
il, for high M, is 
are obtained with 



such procedure being shown in Fig. 

M inferred from ref. j; 
These results for LUp 
f = 0.55eV. 

The steep drop in the resistivity below Tc can be un- 
derstood in a consistent way: the resistivity is dominated 
by spin and phonon scattering at high temperatures until 
the magnetic transition is reached. Below this point the 
negative mobility gap increases along with M, delocal- 
izing a considerable amount of the previously localized 
states. The scaling theory of localization 0| prescribes 
that cr(M) ~ (— A(M))'^. Replacing here the results for 
A(M) of Fig. HI we find that p{Tc)/plT = 0) ~ 20-50, 
in agreement with the values obtained in the experiments 
0, Q|. The negative magnetoresistance (MR) comes 
about as a natural consequence of our microscopic mech- 
anism and further supports our claim that the depen- 
dence of the mobility gap on the magnetization is the 
most relevant factor driving the physics of this material 
near Tc. 

Doped EuBg. Mean-field ana lvsis of Hamiltonian (jSJ 



predicts that Tc scales with (H) 
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. Since cu?, follows 



the same scaling, as stated above, we expect that in the 
series Eui-j^Ca^^Bg the squared plasma frequency should 
scale approximately with Tc, a prediction that would be 
interesting to investigate experimentally with additional 
infrared reflectivity experiments. 

Band structure calculations seem to agree that the con- 
duction band has a strong 5d Eu component. Ca doping 
not only dilutes the magnetic system but also the con- 
ducting lattice. The hoping parameter Uj in eq. ^ is 
then replaced by ti,jPiPj, where = 1 if the site i is oc- 



cupied by a Eu atom and pi = 0, otherwise. The micro- 
scopic problem thus becomes a DE problem in a percolat- 
ing lattice which, at T = K, reduces to a quantum per- 
colation problem Since Ca and Eu are isovalent in 
hexaborides one does not expect the number of carriers to 
depend on x. Nevertheless, since carriers are presumed to 
arise from defects it is difficult to be specific on this issue. 
The mobility edge, on the other hand, is very sensitive 
to the Eu Ca substitution, and should drift toward tjie 
band center. In the paramagnetic regime (T > Tc), Ec 
should move past the Fermi energy at some critical dop- 
ing x^jj, after which the mobility gap becomes positive. 
This determines a crossover from the metallic regime to 
an insulatingbehavior for T > Tc, as seen in the doped 
compounds 0, . At finite T the mobile carriers arise 
from thermal activation across the mobility gap. The re- 
sistivity should display a semiconducting behavior with 
T and its dependence on M should be dominated by an 
exponential factor p(M) ^ exp(A(M)/(fcBT)) Q. Us- 
ing A(M) w Ao(l — aM), as happens in the non-diluted 
case for either of the A(M) ^ situations, we find that 
p{M) ~ exp(— /9M) {f3 is a constant), as seen in the ex- 
periments 
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Figure 3: Schematic phase diagram for Eui-^iCaaiBe in the 
T-x plane in standard notation: PM stands for paramag- 
netic metal, PI for paramagnetic insulator, FM means ferro- 
magnetic metal and FI ferromagnetic insulator. 



As T is lowered below Tc and M increases on the per- 
colating cluster, one expects the crossover from metal- 
lic to semiconducting behavior to occur at larger values 
of x; this is illustrated in Fig. El by the curved dashed 
line separating the FM and FI regions. In the vicinity 
of this line, a sharp metal-semiconducting distinction is 
not possible, resulting possibly in a bad-metal behavior. 
At T = 0, there is a metal-insulator transition occur- 
ring at a concentration xmi > x^jj, which corresponds 
to the quantum percolation transition for a small num- 
ber of carriers (this is different from the usual quantum 
percolation point pq < Pc, which is defined by the lo- 
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calization of all states in the band |4Q|). Even though 
we expect xmi to be close to x^j, the possibility of a 
semiconducting behavior crossing over to metallic at low 
T (for some x^jj < x < xmi) cannot be excluded. Fer- 
romagnetism induced by the DE mechanism is expected 
to persist past xmi as long as the localization length is 
greater than the lattice spacing. Naturally, the critical 
concentration, Xc, where Tc 0, should not be higher 
than pc « 0.69, the classical site-percolation threshold 
for the simple cubic lattice The values of x^jj and 
Xmi in Fig. El vary with carrier density and are expected 
to be sample dependent, since carriers seem to originate 
from defects. In fact, annealing experiments can be quite 
important for the study of the phase diagram. 

Some questions remain as yet unanswered. The value 
of t = 0.55 eV used in FigI211eads to a mean-field estimate 
of Tc w 80 K which is much higher than the exper- 
imental value. At such low carrier concentrations, mag- 
netic polaron formation [l^, must certainly be taken 
into account close to Tc and could lead to a reduction of 
the mean-field estimate. Indeed, Raman scattering stud- 
ies |4^ li^ l close Tc have been interpreted in terms of 
magnetic polarons. A large, positive, MR is observed at 
high fields and T <Tc in EuBg [lEj- H ^ > 0, it is pos- 
sible that high magnetic fields may bring into play states 
of parallel electronic and Eu spins, that were projected 
out to obtain the Hamiltonian Q . One would then have 
to deal with two types of carriers and therefore a positive 
MR. Finally, our model allows for the unusual situation 
of a competition between conduction by thermal activa- 
tion across a small, magnetization dependent, mobility 
gap, and conduction by hopping among localized states. 
Some of the details of the transport properties, such as 
the remarkable scaling of the Hall coefficient with magne- 
tization in Euo.6Cao.4B6 will require a more thorough 
understanding of transport in these circumstances. 

In conclusion, we have proposed a theory for Eu hex- 
aborides, based on the close proximity of the Fermi level 
and a magnetization dependent mobility edge Ec{M). 
This theory arises quite naturally from the Kondo lattice 
problem in the limit of very small number of carriers; it 
accounts, in a quantitative way, for many of the observed 
experimental properties in these materials. Furthermore, 
we also predict a linear scaling of ujp with Tc that can 
be tested in infrared refiectivity experiments. 
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